Dynamics, Rectification, and Fractionation for Colloids on Flashing Substrates 
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We show that a rich variety of dynamic phases can be realized for mono- and bidisperse mixtures 
of interacting colloids under the influence of a symmetric flashing periodic substrate. With the 
addition of dc or ac drives, phase locking, jamming, and new types of ratchet effects occur. In some 
regimes we find that the addition of a non-ratcheting species increases the velocity of the ratcheting 
particles. We show that these effects occur due to the collective interactions of the colloids. 
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The motion and ordering of colloidal particles on two- 
dimensional (2D) periodic substrates has been attract- 
ing growing interest due to recent experimental break- 
throughs that permit the creation of static ID and 2D 
periodic substrate arrays using optical and holographic 
techniques 0, 00 HSU 00 • K is also possible to cre- 
ate dynamic periodic arrays, such as flashing or shifting 
traps 0, 0, . Colloidal particles interacting with 
periodic substrates are ideal model systems for studying 
general problems in condensed matter physics, such as 
the ordering and melting of commensurate and incom- 
mensurate elastic lattices on periodic surfaces. Problems 
of this type can arise in vortex crystal ordering in su- 
perconductors 0| or Bose-Einstein condensates (BECs) 
with periodic pinn ing sites [lij as well as molecules ad- 
sorbed on surfaces |l5j . The dynamics of colloids moving 
over periodic substrates is also relevant to understanding 
depinning phenomena Il6| and phase locking . In 
addition to the scientific interest in these systems, there 
are technical applications for colloids moving over peri- 
odic arrays, such as the fractionation or segregation of 
colloidal mixtures 0, 0, 0]- F° r example, species 
fractionation can be achieved by flowing colloidal mix- 
tures over an array at an angle such that the motion 
of one colloid species locks to a symmetry direction of 
the periodic array while the other species moves in the 
driving direction 0, . Other segregation phenomena 
occur when the substrate is dynamic, and rectification or 
ratchet devices can be constructed in which one species 
ratchets at a different velocity than the other [Io|. New 
ratchet-based logic devices have also been realized 
which may be useful for understanding similar solid state 
nanoscale devices [lij . 

Assemblies of interacting particles have been exten- 
sively studied for systems with random disorder and pe- 
riodic substrates. Collectively interacting particles on a 
flashing substrate is a new class of problem that can be 
realized experimentally using dynamic substrates. For 
particles interacting with optical traps, flashing can be 
achieved simply by modulating the laser power. It should 
be possible to create similar flashing potentials for vor- 



tices in BECs or ions interacting with optical arrays [14( . 

In this work we consider mono- and bidisperse assem- 
blies of charged colloids interacting with flashing 2D sym- 
metric periodic substrates. When a dc drive is applied 
to this system, phase locking in the velocity force curves 
occurs. If an external ac drive is applied instead of a dc 
drive, we find that new types of ratchet effects can be 
realized. For bidisperse colloidal assemblies, the relative 
velocity between the two colloid species is affected by the 
density, number ratio, and charge disparity of the species. 
We find that collective effects play an important role in 
the ratcheting behavior. For a system containing a non- 
ratcheting species A, the introduction of a second species 
B that does ratchet can induce ratcheting in species A 
due to the interactions between the colloid species. The 
effectiveness of the induced ratcheting goes through a 
peak as the fraction of ratcheting particles is changed 
at fixed colloid density. At high densities, the colloids 
reach a jammed state and both species ratchet at the 
same velocity. We also find a remarkable phenomenon for 
strongly interacting particles, where it is possible to in- 
crease the drift velocity of a ratcheting species by adding 
non-ratcheting colloids to the system. This effect oc- 
curs when the disorder introduced by the non-ratcheting 
species breaks apart structures formed by the ratcheting 
species and allows the latter to couple more effectively to 
the substrate. Our results can be tested experimentally 
for colloids interacting with flashing optical traps. 

We consider a 2D Brownian dynamics (BD) simulation 
of N interacting colloids with periodic boundary condi- 
tions in both the x and y-directions. We neglect hydrody- 
namic and excluded volume interactions, which are rea- 
sonable assumptions for strongly charge-stabilized parti- 
cles in the low volume fraction limit. The overdamped 
equation of motion 0] for colloid i is given by 

dTt 

where the damping constant rj is set to un it^- The colloid- 
colloid interaction force is F£ c = q lt J2i^j ~ ^i^( r y Oj 
where the colloid-colloid potential is a Yukawa or 
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1: (a) {V x ) vs F dc for a system with iV/ATp = 1 and 
- 361 for A = 1 at varied flashing period. From bottom 
to top, T s = 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 
600, 650, 700, 750, and 800 BD steps, (b) Same as in (a), 
with a fixed T s = 300 and increasing substrate strength A. 
From top to bottom, A = 1,2, 3, 4, 5, 6, 7, 8, 9, and 10. All 
curves have been vertically offset for clarity. 

screened Coulomb interaction with the form V(rij) = 
Eoiqj/lTi-TjDexpi-Klvi-Tjl). Here E = Z* 2 /(4iree ) 
where charge is measured in units of Z* and e is the sol- 
vent dielectric constant. Lengths are measured in units 
of do, assumed of order a micron, while time and force 
are measured in units of r = rj/E and F = E /a,Q. qju\ 
is the charge on particle 1/n is the screening length 
which is set to 2ao, and ^i(j) is the position of particle 
For monodisperse assemblies, qi/qj = 1, while for 
bidisperse assemblies, qi/qj ^ 1. The thermal force F T 
is modeled as random Langevin kicks with the proper- 
ties (Ff ) = and (F T (i)F T (i')) = 2nk B T8{t - t'). The 
substrate force Ff — Asm.(2Trx/a x )x + B sm{2-Ky / a y )y 
with N p minima is flashed on and off with a square wave 
that has a 50% duty cycle beginning in the off state 
at t = 0. Here a x = 4ao is the substrate periodicity, 
a y = \/3a x /2, A = B is the strength of the substrate, 
T s is the flashing period, and uo s — 2w/T s . For most of 
this work, the sample size is 76ao x 65.8eto. The term 
pext re p resen ts an externally applied drive. For a dc 
drive, Ff xt — F dc ic, and for an ac drive of frequency u>d, 
phase offset <j>, and no force off set from zero, Y\ xt = 
F ac sin((wd + (f>)t)x. We allow the system to reach a 
steady state, which takes less than 5000 BD steps, before 
measuring the time-averaged particle velocities of each 

species, {V A{B) ) = (T av N A ( B ))~ l J2?Ji (B) El=o * ' V *M' 
where r av > T s , r av = 5000 to 1 x 10 6 BD steps. 

We first consider the case of monodisperse colloidal 
assemblies interacting with a flashing substrate and a 
dc drive, where we find phase locking. For a parti- 
cle moving over a periodic substrate with an additional 
ac drive, phase locking takes the form of Shapiro steps 
(lif. in which the particle velocity is periodically modu- 
lated by the dc motion over the substrate at a frequency 
u m = a x rj/(F dc T). When ui m — nujd, with n an integer, 
a synchronization of frequencies occurs and the particle 
motion stays locked at the same dc velocity over a range 
of dc drive, producing steps in the velocity-force curve. 
In our system we consider a dc drive and a flashing sub- 
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FIG. 2: (a) (Vx) vs F ac for a system with uj s /LOd = 1 and 
N p = 64. From top to bottom maximum, N/N v = 1, 1.25, 
1.28, 1.33, 1.375, 1.41, 1.44, 1.47, 1.5, and 2.0. Inset: (V x ) for 
the same system in (a) for varied phase difference <f> between 
the ac drive and flashing substrate, (b) Same as in (a) for a 
higher N p = 361 with N/N p = 1 (top peak), 1.025, 1.05, 1.10, 
1.39, and 1.66 (bottom smooth curve). Inset: the maximum 
value of (V x ) vs N/N p for the system in (b). 

strate, and there is no external ac driving force. The 
synchronization occurs as ui m = nu> s . Here, the flashing 
substrate permits the realization of phase locking with- 
out an external ac drive. 

In Fig. 1(a) we show a series of velocity- force curves 
for varied uj s for the case iV = N p . A depinning thresh- 
old appears at F dc = 0.34 followed by a series of phase 
locking steps. For higher uj s , steps at higher F dc can be 
resolved. In Fig. 1(b) we plot the velocity force curves 
for varied substrate strength A at fixed T s — 300. Larger 
numbers of steps that have a saturated width appear for 
larger A. These saturated steps are distinct from Shapiro 
steps, where the width of the nth step varies with ac am- 
plitude as the nth-order Bessel function. For F dc > A, 
particles are no longer trapped in the minima when the 
substrate cycles on, and much weaker Shapiro steps ap- 
pear. When the number of colloids is incommensurate 
with the substrate, the steps become smooth and have 
a finite slope as additional soliton type pulses of motion 
are present which do not lock to u> s . 

If an ac drive is applied instead of a dc drive, it is 
possible to produce a ratchet effect. Consider the case 
where F ac < A, such that if the substrate potential were 
static, the particle would remain trapped within a single 
pinning plaquette. The particle is able to move under the 
influence of the ac drive during the period of time when 
the substrate flashes off, but will be trapped in whichever 
plaquette it reached when the substrate flashes back to 
the on state. Depending on the amplitude, frequency uj<i, 
and phase of the ac drive relative to the flashing of the 
substrate, the particle may be able to move forward a dis- 
tance na in the x direction with each cycle of the pinning, 
where n is an integer, but is trapped inside a pin during 
part or all the period of time when it would otherwise 
have moved a distance —na, resulting in a net dc motion. 
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FIG. 3: Ratcheting effect for colloidal mixtures in a sam- 
ple with N p = 722. Open squares: (V^) for non-ratcheting 
species A; filled diamonds: (V a f) for ratcheting species B; 
plus symbols: relative velocity, (a) Fixed Na/N =1/2 for 
increasing TV. (b) Fixed N A /N = 1/2, TV = 200, and q B = 3 
for increasing g.4- (c) Fixed TV — 300 and increasing Na/N. 
(d) Fixed TV = 500 and increasing Na /TV. 

In Fig. 2(a) we show the average velocity in the x direc- 
tion (14) vs F ac for a system with w s /cjd = 1, TV p = 64, 
and <fi = 0. For the commensurate case of TV = 64, the 
effective ratchet velocity goes through a series of peaks of 
decreasing height, with the maximum peak occurring at 
F ac = 1.25. Between the peaks, (14) « 0. The de- 
creasing envelope of peak heights with increasing F ac 
arises since the particles are able to depin during part 
of the cycle when the pins have flashed to the on state, 
decreasing the net rectification. At incommensurate fill- 
ings TV 7^ TV p , the colloid positions are more disordered 
due to the colloid-colloid interaction forces and the rec- 
tification peaks become smeared. In Fig. 2(b) we plot 
results from a larger system of 19 x 19 pins at the same 
a x . At N/Np = 1, peaks again appear as a function of 
F ac , while at the higher colloid densities the peaks are 
smoothed. The peak height is constant for N/N p < 1 
since at submatching densities, particle-particle interac- 
tions do not play an important role and the rectification 
that occurs is a strictly single-particle process. In the in- 
set of Fig. 2(b) we show the monotonic decrease for the 
maximum peak value of (V x ) for increasing TV/TV p , which 
can be fitted well to an exponential decay. The ratchet 
effect also depends on the phase between the flashing 
pinning and ac driving force. In the inset of Fig. 2(a) we 
show V x vs <fi for a fixed F ac = 1.0. V x is positive for 
4> = 0, zero near <f> = n/2, and negative for <p = n. 

We next consider mixtures of two colloidal species with 
different charge. The coupling to the ac drive is weaker 
for the species with smaller charge. We first show how 
the density of the system affects the ratcheting for an ac 
drive at which species B with charge qs — 3 ratchets and 
species A with charge q& = 1 docs not. In Fig. 3(a) we 
plot the average velocities (V^) and (Vj 3 ) vs N, with 
fixed N A /N = 1/2, A = 1.0, T s = T d = 500, and 
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FIG. 4: Snapshot of colloid positions (circles) and trajectories 
(lines) in a portion of the sample for a fixed period of time 
for the system in Fig. 3(a) at TV = 40. Dark circles are the 
ratcheting species B and light circles are the non-ratcheting 
species A. (b) Snapshot of the jammed system in Fig. 3(a) 
at TV = 700. (c) Snapshot at N A /N = 0.2 in Fig. 3(d). (d) 
Snapshot at N A /N = 0.6 in Fig. 3(d). 

Qb/qa = 3. For low densities, species B moves at a finite 
velocity while (V^) is close to zero. As TV increases, the 
particles interact more strongly and the ratcheting par- 
ticles B push the A particles forward, inducing a ratch- 
eting effect for species A. In Fig. 4(a) we illustrate the 
motion of the colloidal mixtures at a low density. All of 
the species B particles are mobile, as indicated by the 
trajectories. In several places species A particles can be 
seen to be immobile. However, in the middle and right 
portions of the figure, two A particles have trails indi- 
cating motion that was produced as a species B particle 
pushed the species A particle. This is the mechanism 
which causes the ratcheting in species A. We note that 
if only species A is present, no ratchet effect occurs for 
any value of TV unless F ac is increased. The removal of 
the last few B particles slows the system in a non-linear 
fashion until the whole system stops at Na/N = 1. 

As the total number of particles increases for fixed 
Na/N = 1/2, shown in Fig. 3(a), (V^) drops as more of 
the ratcheting energy is transfered to species A. We find 
a maximum m (V/) at TV w 400. The velocity of both 
species decreases for TV > 400. For TV > 700 the system 
jams and the particles can no longer exchange positions, 
so both species move with the same velocity. In Fig. 4(b) 
at N = 700, the species tend to clump together, with 
species B tending to form crystallites. The jamming or 
freezing transition we find is similar to the freezing at 
increased density observed in systems of two species of 
particles moving in different directions [lif. 



(b) 



, o o o oo o • „o o • o • o 

v • • • ; * • °°i 



(d) 



4 



In Fig. 3(a), the relative velocity {V^} — {V^} between 
species A and B is a measure of the fractionation, and it 
drops to zero when the system jams. In Fig. 3(b) we plot 
(V A ) and (V^) vs varied charge qA for the same system 
in Fig. 3(a) at fixed density with N = 200. Here (V X A ) 
increases monotonically with qA, and for qA > 2.5, the 
velocities of the two species are the same. (V^) decreases 
with increasing qA for qA < 1.0 since species A needs 
to be pushed in order to ratchet for these values of qA, 
and the coupling between the species increases with qA, 
resulting in a larger drag on species B. For qA > 1, the 
A species starts to ratchet on its own and {V^} recovers 
as the drag from species A diminishes. 

In Fig. 3(c) for a low density at N = 300 we show 
V x as a function of Na/N. As the fraction of ratcheting 
colloids goes to zero, the velocity of both species is mono- 
tonically reduced. In the low density regime the system 
is always liquidlike. In Fig. 3(d) we show a very different 
behavior at a much higher N = 500. Here, increasing the 
fraction of non-ratcheting particles Na/N increases the 
average velocity of both the ratcheting and non-ratchcting 
species. The relative velocity between the two species 
also increases for higher Na/N. These effects are more 
pronounced for higher N. For small Na/N at high den- 
sities, the colloid-colloid interactions dominate and the 
system acts as a rigid unit so that all the colloids must 
ratchet at the same velocity together, as in Fig. 4(c). 
In general some colloids fall at locations that are incom- 
mensurate with the substrate, so they do not ratchet as 
effectively and slow the entire system. In Fig. 3(d) for 
low Na/N, (V a ) = (Vj 3 ), indicating that the system is 
moving elastically. For higher fractions the velocities dif- 
fer, indicating that tearing of the lattice occurs, and at 
the same time the velocities of both species increase. We 
observe that in this case the system is much more disor- 
dered and the colloids act more as single particles rather 
than as a collective unit. Thus, the addition of the non- 
ratcheting particles effectively breaks up the system or 
adds disorder, as in Fig. 4(d), which allows the system to 
behave as a liquid rather than a solid. This improves the 
coupling of the ratcheting colloids to the substrate since 
they are no longer held in incommensurate positions, and 
increases the overall effectiveness of the ratcheting. We 
have also performed extensive simulations for other pa- 
rameters such as varied ratios of A and B, changing the 
steepness of the substrate potentials, densities, Na/N, 
and T. In all cases we find the same general features of 
the ratchet effects described here, leading us to believe 
that the effects are very generic. 

To summarize, we have shown that a variety of dynam- 
ical behaviors, including phase locking, ratchet effects, 
and jamming, are possible for colloids and colloidal mix- 
tures interacting with flashing symmetric periodic sub- 
strates. For a flashing substrate with a dc drive, a phase 
locking occurs that is distinct from the typically studied 
Shapiro step phase locking. If a strictly ac drive is applied 



to a flashing substrate, a ratchet effect can occur, with a 
ratcheting direction that depends on the relative phase of 
the two ac signals. For mixtures of particles in regimes 
where one species ratchets and the other does not, we 
show that the ratcheting species can induce a ratchet ef- 
fect in the non-ratcheting species, at the cost of a reduced 
effectiveness of ratcheting in the first species. For high 
densities the system jams and both species move at the 
same velocity. We have demonstrated the ratchet effect 
for a wide range of parameters for the mixtures. In the 
dense regime we find that adding nonratchcting colloids 
can increase the drift velocity of the ratcheting species. 
This is due to the non-ratchcting particles breaking up 
the elastic flow of the ratcheting species and improving 
the coupling to the substrate. Our results should be read- 
ily testable for colloids on dynamical optical traps. 
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